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Implications from the optical to UV flux ratio of Fell 
emission in quasars 
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ABSTRACT 

We investigate Fen emission in Broad Line Region (BLR) of AGNs by analyzing 
the Fen(UV), Fcn(A4570) and Mgn emission lines in 884 quasars in the Sloan Dig- 
ital Sky Survey (SDSS) Quasar catalog in a redshift range of 0.727 < z < 0.804. 
Fen(A4570)/Fen(UV) is used to infer the column density of Fe n-emitting clouds and 
explore the excitation mechanism of Fe n emission lines. As suggested before in various 
works, the classical photoionization models fail to account for Fen(A4570)/Fen(UV) 
by a factor of 10, which may suggest anisotropy of UV Fe n emission; otherwise, 
an alternative heating mechanism like shock is working. The column density distri- 
bution derived from Fen(A4570)/Fen(UV) indicates that radiation pressure plays 
an important role in BLR gas dynamics. We find a positive correlation between 
Fen(A4570)/Fen(UV) and the Eddington ratio. We also find that almost all Fen- 
emitting clouds are to be under supcr-Eddington conditions unless ionizing photon 
fraction is much smaller than that previously suggested. Finally we propose a physical 
interpretation of a striking set of correlations between various emission-line properties, 
known as "Eigenvector 1" . 

Key words: galaxies: active - galaxies: nuclei - quasars: emission lines - line: for- 
mation - atomic processes - radiation mechanisms: general 



1 INTRODUCTION 

According to the models of nucleosynthesis, much of Fe 
comes from Type la supernovae (SNe la), while a elements 
such as O and Mg come from Type II supernovae (SNe II). 
Because of the difference in lifetime of the progenitors, Fe- 
enrichment delays relative to a elements. Hence the abun- 
dance ratio of Fe to a elements [Fe/a] should have a sud- 
den break at 1-2 Gyr after the initial burst of star forma- 
tion (Hamann & Ferland 1993; Yoshii, Tsujimoto & Nomoto 
1996; Yoshii, Tsujimoto & Kawara 1998, but see also re- 
cent studies, such as Matteucci et al. 2006 and Totani et 
al. 2008, indicating a significant number of SNe la on rel- 
atively short timescales). Under the assumption that the 
Fell/Mgll flux ratio reflects [Fe/Mg], various groups have 
measured Fell/Mgll in high-redshift quasars hoping to dis- 
cover such a break (e.g., Elston, Thompson & Hill 1994; 
Kawara et al. 1996; Dietrich et al. 2002, 2003; Iwamuro et 
al. 2002, 2004; Freudling, Corbin & Korista 2003; Maiolino 
et al. 2003; Kurk et al. 2007; Sameshima et al. 2009). How- 
ever, these efforts have ended up with finding a large scatter 
of Fell/Mgll showing little evolution. 
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A doubt is, thus, cast on the assumption that Fe Il/Mg II 
reflects [Fe/Mg]. For examples, Verner et al. (2003) sug- 
gested that Fell/Mgll depends not only on the abundance 
but also on the microturbulence of Fe II-emitting clouds. 
Tsuzuki et al. (2006) showed that Fell/Mgll correlates with 
the X-ray photon index, the full width at half maximum 
(FWHM) of Mgll, the black hole mass, etc. For these rea- 
sons, prior to deriving the abundance from Fell/Mgll, we 
must first clarify these non-abundance effects on the Fell 
emission as well as the source of the Fe II excitation. 

Column densities of clouds, which are considered to be 
one of the non-abundance factors largely affecting the Fell 
emission, have lately attracted attention for their signifi- 
cances in determining whether or not the radiation pressure 
plays an important role in BLR gas dynamics. Marconi et al. 
(2008) considered the effect of radiation pressure from ion- 
izing photons on estimating of the black hole mass, which 
is based on the application of the virial theorem to broad 
emission lines in AGN spectra, and suggested that the black 
hole mass can be severely underestimated if the effect of 
radiation pressure is ignored. Netzer (2009) then used the 
Mbh — 0"* relation for a test of this suggestion, where Mbh 
is the black hole mass and a, is the velocity dispersion of 
host galaxies, concluding radiation pressure effect is unim- 
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Figure 1. Simplified Fc II Grotrian diagram. Note that each Level 
represents a large number of levels that have nearly the same 
excitation energies. Solid arrow indicates UV Fe II emission lines, 
while dashed arrow indicates optical Fell emission lines. 



portant, while Marconi et al. (2009) found the importance 
of radiation pressure by taking into account the intrinsic 
dispersion associated with the related parameters, in partic- 
ular, column densities of BLR clouds. However, there are no 
reliable column density estimates from observations up to 
date. 

In this paper, we will estimate column densities of 
quasars selected from the SDSS using Fell, for investigating 
the excitation mechanism of Fe II emission and the BLR gas 
dynamics. In section 2, we perform numerical calculations 
of Fe II emission lines to establish a method for estimating 
column densities. In section 3, Fell emission lines in the UV 
and optical as well as Mgll emission lines are measured in 
the SDSS quasars. The results and discussion about Fell 
emission mechanism, radiation pressure, and the variety of 
quasar spectra called as "Eigenvector 1" are given in sec- 
tion 4. Throughout this paper, we assume a cosmology with 
Q m = 0.3, Ov = 0.7 and H = 70 km s" 1 Mpc" 1 . 



2 METHODOLOGY 

In the following, we use Fell(UV) to denote the UV Fell 
emission lines in 2000 < A < 3000 A, FeIl(A4570) to the 
optical Fell emission lines in 4435 < A < 4685 A, and Mgll 
to the Mgll A2798 emission line. 



2.1 Fen(A4570)/Fen(UV) to measure column 
densities 

Figure 1 shows a simplified Fell Grotrian diagram. As can 
be seen, Level 3— Level 2 transitions give rise to optical Fell 
emission lines such as the FeIl(A4570) bump. Branching ra- 
tio of Level 3— Level 1 UV resonance transitions is signifi- 
cantly higher than that of Level 3— Level 2 optical transi- 
tions. Hence strong optical Fell emission requires a large 
optical depth between Level 1 and Level 3 such as T13 > 10 3 
in order to transform UV Fell lines to optical Fell lines 
through a large number of scatterings (cf. Collin & Joly 
2000). Thus FeIl(A4570)/Fen(UV) flux ratio must strongly 



depend on T13 and T23 , which are the optical depths for pho- 
tons emitted through Level 3— Level 1 and Level 3— Level 
2 transitions, respectively. If so, the Fe Il(A4570)/Fe Il(UV) 
can be an indicator of the column density. 

Here we use a quite simple model to indicate the de- 
pendence of FeIl(A4570)/FeIl(UV) on the column density. 
First, we ignore the Level 4 in Figure 1 and consider the 
Fell as three-level system. Second, we assume thermal equi- 
librium population between Level 1 and Level 2. Third, we 
assume an expression of the local escape probability given 
by Netzer & Wills (1983) as e tJ = (l-Tij)/nj. Although the 
model adopting these assumptions is obviously too oversim- 
plified, it is useful to qualitatively understand how the line 
ratio depends on the physical parameters. The flux ratio is 
then written as 



Fen(A4570) 
Fell(UV) 



713^31 £31/1^13 

'2.9 eV\ 1 



exp 



/ 2.9 eV \ 
V kT ) 



1 



(1) 
(2) 



where 713 is the population of level 3, and Aij is the sponta- 
neous emission rate from the level i to the level j. Following 
the fact that T13 ^> 1, the ratio can be roughly reduced to 
FeIl(A4570)/FeIl(UV) tx l-e _T23 . Since r 23 is proportional 
to the column density, the ratio is an increasing function of 
the column density. We will discuss further on this matter 
in the following. 



2.2 Photoionization models in the LOC scenario 

Here we will show more sophisticated model calculations 
of FeIl(A4570)/FeIl(UV). We performed photoionization 
model calculations with version C06.02 of the spectral sim- 
ulation code Cloudy, last described by Ferland et al. (1998), 
combined with a 371 level Fe + model (up to ~11.6 eV, 
Verner et al. 1999). The incident continuum is defined as 



/„ oc v auv exp( — hv/kTsB) exp(—kTm/hi>) + av a 



(■■>>) 



The first term in the right-hand side of the equation (3) 
expresses an accretion disk component, which is usually 
called Big Bump. The KTbb and the kTm indicate the 
higher and the lower cut off energies, respectively. The sec- 
ond term expresses a power-law X-ray component, which is 
set to zero below 1.36 eV, while to fall off as v~ z above 
100 keV. The coefficient a is set to produce the optical 
to X-ray spectral index a ox ■ We set these parameters to 
(a uv , a x , a ox , T BB , kT IR ) = (-0.2, -1.8, -1.4, 1.5 x 
10 5 K, 0.136 eV), which are adopted in Tsuzuki et al. (2006). 
This incident continuum illuminates a single cloud with hy- 
drogen density uh, ionization parameter U(= Q/uhc, where 
$ is the ionizing photon flux, c is the velocity of light), 
column density Nh and solar abundance. We calculated 
the models in a range of Nh = W 2 



n H = 



10' 



10 14 cm J and U = 10" 



10 u 



In locally optimally emitting clouds (LOCs; Baldwin et 
al. 1995), each line is emitted from clouds with a wide range 
of gas hydrogen density and distance from the central con- 
tinuum source, and the observed spectra is reproduced by 

1 The optical to X-ray spectral index a ox is defined as /i/(2 
kcV)//„(2500A) = 403.3 Q °- 
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integrating these clouds with an appropriate covering frac- 
tion distribution. The observed emission line flux is, then, 
expressed as: 



0.100 



Llirie. 



F{r, n H ) j '(r)g(n H )dn H dr 



(4) 



where F(r, uh) is the emission line flux of a single cloud at 
a distance r from the central continuum source and with 
uh, f(r) is a cloud covering fraction with distance r, and 
g(riH) is a fraction of clouds with nn- Matsuoka et al. (2007) 
showed that O I and Call emission lines in quasars, which are 
likely to emerge from the same gas as the Fe II emission lines, 
are well reproduced by a LOC model with /(r) oc r _1 and 
g(riH) oc nj^ , hence we have adopted these covering distri- 
butions. Baldwin et al. (1995) suggested that since clouds at 
large distances from the continuum source will form graphite 
grains which heavily suppress the line emissivity, clouds with 
$ < 10 18 s -1 cm -2 must be excluded from integration. Thus 
we have calculated equation (4) for clouds which correspond 
$ > 10 18 s^cm" 2 . 

Figure 2 illustrates Fe Il(A4570)/Fe Il(UV) as a func- 
tion of the column density for photoionized clouds. 
Fe Il(A4570)/Fe Il(UV) increases as the column density in- 
creases. It also shows that FeIl(A4570)/FeIl(UV) decreases 
as the microturbulent velocity increases, consistent with 
Verner et al. (2003), and that FeIl(A4570)/FeIl(UV) is max- 
imum when no microturbulence is assumed to exist. This 
is explained as follows. Increasing microturbulent velocities 
broadens the line absorption profile, resulting in enhance- 
ment of the continuum photoexcitation. This effect is rela- 
tively large for high energy levels where the collisional exci- 
tation is inefficient for their high excitation potential. Thus 
large microturbulent velocities relatively enhances the con- 
tinuum photoexcitation to the Level 4 in Figure 1, lead- 
ing to emit more fluxes in the UV Fell emission lines, thus 
decreasing Fe Il(A4570)/Fe Il(UV). For photoionized clouds, 
Fe Il(A4570) /Fe Il(UV) can be used as a column density indi- 
cator unless microturbulent velocities vary much from cloud 
to cloud. 

It is noted that calculations adopting the other 
shapes of the incident continuum 2 show little changes for 
FeIl(A4570)/FeIl(UV), indicating little dependence on the 
spectral energy distribution (SED) of ionizing continuum. 



2.3 Collisionally ionized models 

As an alternative to photoionization models, we consider 
models in which a cloud is assumed to be in collisional equi- 
librium at a given electron temperature T e , and call them 
collisionally ionized models. In these models, the specific 
heating mechanism is not accounted and arbitrary electron 
temperatures are given. It is noted that Fe II is collisionally 
excited in either mechanically heated clouds such as through 
shocks or photoionized clouds. In the latter case, the heating 
mechanism is specified to be the eating of the incident UV 
and X-ray photons. 



2 Two types of SED given by Nagao, Murayama & Taniguchi 
(2001) arc adopted. One is set to reproduce the ordinary SED of 
broad-line Seyfert 1 galaxies, and the other is set to reproduce 
that of narrow-line Seyfert 1 galaxies. 
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Figure 2. Fe II(A4570)/Fc II(UV) vs. column densities based on 
the photoionization models. V tur b is the microturbulence veloc- 
ity of clouds. Note that, in the LOC, the calculation results are 
not dependent on njj and U since clouds with different hydro- 
gen density and distance from central continuum source are all 
integrated. 



Joly (1987) offers collisionally ionized model calcula- 
tions. In her models, Fell is approximated by a 14-level (up 
to ~5.7 eV), and the emission region is assumed to be a ho- 
mogeneous slab with constant hydrogen density nH , column 
density Nh, electron temperature T e , and not to receive any 
external radiation. The electron temperature ranges from 
6,000 K to 15,000 K, the density from 10 10 to 10 12 cm" 3 
and the column density from 10 21 to 10 25 cm -2 . 

Figure 3 shows FeIl(A4570)/FeIl(UV) as a func- 
tion of the column density, taken from Joly (1987). 
As expected from equation (2), it can be seen that 
FeIl(A4570)/FeIl(UV) increases with the column density 
while decreases with temperature. However, Figure 3 shows 
that FeIl(A4570)/FeIl(UV) does not so much depend on 
the temperature in 6,000 < T e < 10,000 K. It is noted 
that Collin et al. (1980) and Joly (1987) showed emission 
line ratios including Fell in quasars are well accounted by 
cold clouds with 6, 000 < T e < 10, 000 K. Thus it is rea- 
sonable to assume that Fe II-emitting clouds have the tem- 
peratures in a range of 6, 000 < T e < 10, 000 K, and the 
FeIl(A4570)/FeIl(UV) depends little on the temperature 
while strongly varies with the column density. We fit the 
data of 6, 000 < T e < 10, 000 K models by a 4th-order poly- 
nomial, and find 

y = -0.72 + 0.28x - 0.046x 2 + 0.031x 3 - 0.0089a; 4 (5) 



= log Feii(A4570)/Fen(UV) and x = 
2 1 — 23. We will use this relation to estimate 



where y 
log Nh [cm 
the column density of quasars 



3 ANALYSIS OF QUASAR SPECTRA 
3.1 Sample selection 

We have analyzed quasar spectra selected from the fourth 
edition of the SDSS Quasar catalog (Schneider et al. 2007). 
SDSS uses a dedicated 2.5 m telescope at the Apache Point 
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Figure 3. Fc II(A4570)/Fc II(UV) vs. column densities based on 
the collisionally ionized models taken from Joly (1987). The 
density ranges from 10 10 to 10 12 cm -3 and the column den- 
sity from 10 21 to 10 25 cm . Circles indicate the models with 
6, 000 < T e < 10, 000 K, triangles with T e = 12, 500 K and di- 
amonds with T e = 15, 000 K. The solid line indicates the fitted 
4th-ordcr polynomial for samples with 6, 000 < T e < 10, 000 K. 
The dotted and the dashed lines connect the median values of 
Fe II(A4570)/Fe II(UV) at each column density for samples with 
T e = 12, 500 K and T e = 15, 000 K, respectively. 



Observatory equipped with a CCD camera to image the sky 
in five optical bands, and two digital spectrographs, one cov- 
ering a wavelength range of 3800 A to 6 150 A and the other 
from 5800A to 9200A. The spectral resolution ranges from 
1850 to 2200. The fourth edition of the SDSS Quasar cat- 
alog consists of the objects in the Fifth Data Release, and 
contains 77,429 quasars. 

In order to measure the Fell(UV) and FeIl(A4570) 
emission lines simultaneously, we have selected spectra 
which cover the wavelength range of 2200 to 5100 A in the 
rest frame, corresponding to the redshift range from 0.727 
to 0.804. 2,189 objects meet this requirement. Then we have 
checked the signal to noise ratio (S/N) per pixel for each 
spectrum which fulfills median S/N > 10 per pixel at their 
continuum levels for accurate flux measurements. 946 ob- 
jects meet this requirement. All the spectra were inspected 
by eyes and 62 spectra were rejected because of wavelength 
discontinuity, terrible contamination by host galaxy star 
light, etc. Our final sample thus consists of 884 spectra. 

Prior to the measurements, the quasar spectra were 
dereddened for the Galactic extinction according to the dust 
map by Schlegel, Finkbeiner & Davis (1998) using the Milky 
Way extinction curve by Pei (1992). 

Since SDSS quasars are flux-limited in the survey, low 
luminosity quasars are lost in the SDSS sample. Figure 4 
shows luminosity versus redshift for our sample. Because of 
the narrow redshift coverage for our sample, the minimum 
luminosity is almost constant and roughly estimated to be 
ALsioo = lO 44 ' 7 ergs/s. 



3.2 Continuum and line fitting 

In the UV to optical, the quasar continuum is composed 
of (i) the power-law continuum Fjf , (ii) the Balmer con- 
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Figure 4. Luminosity AL5100 vs. redshift for our sample. Lumi- 
nosity distribution is also displayed in the right panel. The dashed 
line shows a rough estimate of minimum quasar luminosity for our 
sample. 



tinuum F^ aC and (iii) the Fell pseudo-continuum F^ eI1 . 
Thus we assumed a following formula as a model continuum 

pcont . 



A = f\ + r\ + #\ 



3.2.1 Power-Law continuum 

The power-law continuum is simply written as follows: 
A 



(6) 



F\ L = F5100 



15100/ 



(7) 



The free parameters of this model are a scaling factor F5100 
and a power-law index a. We chose three fitting ranges, 
2200-2230 A, 4180-4220 A and 5050-5100 A, as continuum 
windows, since these area have little emission lines (see Fig- 
ure 5). There are, however, the Balmer continuum and the 
Fell pseudo-continuum underneath these regions, requiring 
some corrections. 

Tsuzuki et al. (2006) gives 14 quasar spectra covering a 
wide wavelength range and measured accurately their con- 
tinuum levels. We fitted power-law continuum models to 
their spectra in the continuum windows, and compared the 
continuum levels with those given by Tsuzuki et al. (2006). 
We found that our method systematically overestimates the 
continuum levels, 10.1% at 2200-2230 A, 5.7% at 4180-4220 
A and 3.4% at 5050-5100 A. According to these results, we 
first reduced the flux densities of the object by these amount 
at each continuum window, then fitted the power-law con- 
tinuum. An example of the fitted power-law continuum is 
indicated as the dashed line in Figure 5. The measurement 
error of the continuum levels is estimated to be less than 
10%. 



3.2.2 Balmer continuum 

Grandi (1982) gives a formula describing the Balmer contin- 
uum produced by a uniform temperature, partially optically 
thick cloud: 
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(8) 



SDSS J1 14343.06 + 674455.5 



where B\(T e ) is the Planck function at the electron tem- 
perature T e , and tbe is the optical depth at the Balmer 
edge at A = 3646 A. Kurk et al. (2007) assumed gas clouds 
of uniform temperature (T e = 15, 000 K) and the optical 
depth fixed to tbe = 1, and fit equation (8) to their sam- 
ple quasar spectra to estimate the strength of the Balmer 
continuum (see also Dietrich et al. 2003). We followed their 
method and assumed T e — 15, 000 K, tbe = 1. The only one 
parameter, namely the scale factor FsaC, is set free and is 
decided by fitting equation (8) to the power-law subtracted 
spectrum at 3600—3645 A. An example of the fitted Balmer 
continuum is indicated as the dotted line in Figure 5. 

3.2.3 Fe II pseudo- continuum 

Since Fe II has enormous energy levels, neighboring emission 
lines contaminate heavily with each other, which makes it 
difficult to measure the Fe II emission lines. One approach to 
measure the Fe II emission lines is to use Fe II templates. So 
far, several Fell templates are derived from the narrow-line 
Seyfert 1 galaxy, I Zw 1. 

In the UV, Vestergaard & Wilkes (2001) and Tsuzuki 
et al. (2006) give their Fell templates. The template given 
by Vestergaard & Wilkes (2001) do not cover around Mgll 
line. Tsuzuki et al. (2006) used a synthetic spectrum cal- 
culated with the Cloudy photoionization code in order to 
separate the Fell emission from the Mgll line, and derived 
semiempirically the Fell template which covers around the 
Mg II line. Since we want to measure the Mg II emission line, 
we decided to use the UV Fell template given by Tsuzuki 
et al. (2006). 

In the optical, Veron-Cetty, Joly & Veron (2004) and 
Tsuzuki et al. (2006) open their Fe II templates to the public. 
Veron-Cetty et al. (2004) carefully analyzed the Fell emis- 
sion lines in I Zw 1, finding that the Fell lines are emitted 
from both BLR and Narrow Line Region (NLR). They suc- 
ceeded to separate them and called the broad line system 
LI and the narrow line system N3, respectively. Tsuzuki et 
al. (2006) also analyzed the spectrum of I Zw 1 and derived 
the optical Fe II template, which was however not separated 
into the BLR and the NLR components. We applied both 
the broad line system LI template given by Veron-Cetty et 
al. (2004) and the optical Fell template given by Tsuzuki et 
al. (2006) to all of our samples, finding that the latter has 
a slightly smaller average Xv value (median xt ~ 1-48 for 
Tsuzuki et al. 2006, while median xt ~ 1-58 for Veron-Cetty 
et al. 2004). Here we adopt to use the optical Fell template 
given by Tsuzuki et al. (2006). 

Prior to applying the Fell template to each quasar, 
broadening of the template spectrum is needed. Thus we 
modeled the Fell flux density as follows: 
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Figure 5. Top: Sample quasar spectrum (gray line) with the 
fitted power-law continuum F^ L (dashed line), the Balmer con- 
tinuum F^ aC (dotted line), the Fell pseudo-continuum F^ eI1 
(dot-dashed line) and the sum of the three continua F^° nt (solid 
line). Arrows indicate the continuum windows adopted in the 
power-low continuum fitting. Horizontal thick bars indicate the 
masked region adopted in the Fe II pseudo-continuum fitting. Bot- 
tom: The continuum subtracted spectrum. 
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where x = In A, c is the velocity of light and FWHM co „„ rep- 
resents the FWHM of the convolved Gauss function. We first 
calculated the equation (9) for FWHM co „„ = - 5000 km/s 
stepped by 100 km/s, thus prepared 51 Fell emission line 
models. Then we flux-scaled each model to fit the contin- 
uum subtracted spectrum (i.e., the spectrum after subtract- 
ing the power-law and the Balmer continuum), and adopted 
the model which gives the smallest xt value. We decided 
mask regions as follows; 2280-2380A for C II] A2326; 2400- 
2480A for [Neiv] A2423 and Felll; 2750-2850A for Mgll; 
3647-4000A for high-order Balmer lines; 4050-4150A for RS; 
4300-4400A for H 7 ; 4600-5050A for Hell A4686, H/3 and 
[Oiii] AA4959,5009. An example of the fitted Fell pseudo- 
continuum is indicated as the dot-dashed line in Figure 5. 



3.2.4 Mgll lines 

After subtracting the continuum components, we have mea- 
sured the Mgll emission line for estimating the black hole 
mass from its FWHM. We have fitted the Mgll emission 
line profile by two gaussian components. Figure 6 shows an 
example of the Mg II emission line fitting. 
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Figure 6. Mgll emission line {gray line) with fitted two gaussian 
components {dashed lines) and the sum of the two gaussians {solid 
line) . 



3.3 Black hole mass and Eddington luminosity 

For the classical black hole mass estimate in which the radia- 
tion pressure effect is neglected, we use the following formula 
given by McLure & Jarvis (2002): 

, j . n at r„..„, . ,-\2 

I A-L3000 



M 



3.37 



V 10 37 W 



FWHM(Mgn) 



km s 



M-. 



(10) 



The classical Eddington luminosity is given as follows: 



4ircGMBH,om p 



1.26 x 10" 



; MbH.Q 

Mq 



[ergs/s] (11) 



where c is the velocity of light, m p is the proton mass, and 
<tt is the Thomson cross-section. 

On the other hand, Marconi et al. (2008) suggested that 
the force of the radiation pressure should be corrected to de- 
rive the black hole mass. They give the radiation pressure 
corrected black hole mass MBH,rad and Eddington luminos- 
ity L E dd,rad as follows: 

Lboi 



M B H,rad = MbH.O 



LiEddfl 
LEddfi 



a + 



cftNe 



a + a/{aTN H ) 



Lion 

Lboi 



(12) 
(13) 
(14) 



where L^oi is the bolometric luminosity, Li on is the total 
luminosity of the ionizing continuum (i.e., hv > 13.6 eV), 
and a is the ionizing photon fraction. The second term in 
the right hand side of equation (12) represents the correction 
term of the radiation pressure. We here adopt a bolometric 
correction Lboi = 9AL5100 given by Kaspi et al. (2000). 

3.4 Error estimate 

We performed a Monte-Carlo simulation similar to those 
done in Hu et al. (2008) for estimating the measurement 
errors. The detail of the procedure is as follows. 

(i) Generating a composite spectrum. Following Van- 
den Berk et al. (2001), we generated a composite spectrum 



using all of our samples. This composite spectrum repre- 
sents a typical quasar spectrum for our samples, (ii) Obtain- 
ing typical emission line profiles. We applied the measure- 
ment methods written in §3.2 to the composite spectrum, 
and obtained typical emission line profiles for Fell(UV), 
FeIl(A4570) and Mgll. (iii) Making artificial spectra. We 
combined these line profiles with the power-low continuum 
and the Balmer continuum. Thus the simulated spectrum is 
written as follows. 



A 



F[ L (F 5W0 ,a) + F? aC {F BaC ) 
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[ eIIuv (EW FeIluv ] 

FeII X 457Q 
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+ F 



Mgll 



™{EW FeIIxi5m ) 

(EWMgll, FWHMugll) 



(15) 



Note that, for simplicity, we ignored the broadening of the 
pseudo-Fell continuum. Values of input parameters, which 
are given in the parentheses in equation (15), are randomly 
sampled from probability distributions that are made to re- 
flect the observations. Thus, we generated 1,000 simulated 
spectra, (iv) Generating a noise template. Using all the noise 
spectra produced by the SDSS pipeline for our samples, we 
generated a composite spectrum following Vanden Berk et 
al. (2001) and named it a noise template. This noise tem- 
plate is scaled so that the resulting median S/N to be 10 
per pixel at the continuum level for each simulated spec- 
trum, and is treated as its noise. 

Now we have the 1,000 simulated spectra with their 
noise. The measurement methods written in §3.2 are ap- 
plied to these simulated spectra. We calculate the value 
5 sim = {Pout — Pin) /Pin for each simulated spectrum where 
Pin represents the input parameters (i.e., the values given 
in the parentheses in equation (15)) and Pout represents the 
corresponding measured values for the simulated spectra. 
We regard a S i m , a standard deviation of S S i m , as la error 
of the measurement. Thus we evaluate the measurement er- 
rors to be 16.4% for EW of Fell(UV), 22.9% for EW of 
Fen(A4570) and 7.9% for FWHM(Mgll). The simulation 
implies the measurement error to be 2.9% for the luminos- 
ity, which is less than 10% estimated in the power-law fitting. 
Therefore we decided to evaluate the measurement error to 
be 10% for AL3000 and ALsioo- 



4 RESULTS AND DISCUSSION 

4.1 On the excitation mechanism of Fell emission 

Figure 7 shows the observed Fe Il(A4570)/Fe il(UV) distri- 
bution. As can be seen from the comparison between this 
and Figure 2, our photoionization models underpredict the 
FeIl(A4570)/FeIl(UV) by a factor of 10, failing to account 
for the observations. This result is consistent with the pre- 
ceding study by Baldwin et al. (2004). Additional microtur- 
bulence to the photoionized clouds gets the situation even 
worse. Thus, Figure 7 seems to challenge classical photoion- 
ized pictures of Fe II-emitting clouds. On the other hand, in 
the case of the collisionally ionized models shown in Fig- 
ure 3, the observed FeIl(A4570)/FeIl(UV) flux ratios are 
well reproduced with 10 22 < Nh < 10 24 cm -2 and with 
6, 000 < T e < 10, 000 K. 

These results give two remarks: (1) the Fe II-emitting 
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clouds in quasars are heated to 6, 000 < T e < 10, 000 K; (2) 
the UV and the X-ray photons, which are the heating source 
in our photoionization models, fail to heat the gas to such 
temperatures (probably heat the gas too hot!). One possible 
interpretation is that the Fe II-emitting clouds are heated by 
an alternative mechanism such as through shocks. Here we 
note that there is a reverberation mapping study implying 
shock heating for Fell emission. Kuehn et al. (2008) ana- 
lyzed optical Fell emission bands in the Ark 120, finding 
that they do not respond to the continuum variation. Thus 
the optical Fe II-emitting region may be heated by other 
mechanism than photoionization. These results favor the 
shock heating for the optical Fe II-emitting region, but there 
are also difficulties. First, the amount of shock-processed 
matter would probably be too large. Second, as Kuehn et al. 
(2008) showed, collisionally ionized models failed to match 
the shape of the optical Fe II-emission band. Third, the fact 
that there is no response to the continuum variation for op- 
tical Fe II emission bands can also be interpreted as that the 
emitting region is too large to vary optical Fe II emission in 
observable timescales. Unless shocks are a viable solution, 
the failure of the photoionization model simply indicates 
that it is not predicting the correct heating rate, or that the 
radiative transport calculations are not correct. 

One possible cause disturbing classical photoionization 
models to reproduce the observations may be the assump- 
tion that the Fell emission is isotropic. Ferland et al. (2009) 
recently suggested that UV Fe II lines are beamed toward a 
central source while optical Fell lines are emitted isotropi- 
cally. Then photoionization models can reproduce the ob- 
served UV to optical Fell flux ratio if the Fe II-emitting 
clouds are distributed asymmetrically so that we mainly ob- 
serve their shielded faces. However, this needs special ge- 
ometrical distributions like Type II AGNs; a thick Fe II- 
emitting gas surrounding the central source with a sub- 
stantial covering factor, and intervening between the cen- 
tral region and our eyes. At the present time, it is not clear 
whether or not photoionization models can reproduce the 
emission line strengths other than Fell under such the sit- 
uation. Much broader exploration of photoionization model 
calculations is certainly needed. 
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Figure 7. Observed Fe II(A4570)/Fe II(UV) flux ratio distribu- 
tion. The average of log Fe II(A4570)/Fe II(UV) is -0.8 and the 
standard deviation is 0.2 dcx, in disagreement with the range that 
the photoionization models predict. 
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Figure 8. Estimated column density distribution. The av- 
erage and the standard deviation of the distribution are 
(log Nh , o"i og pf H ) = (22.8,0.5), providing observational support 
for the assumption adopted in Marconi et al. (2009). 



4.2 Column density distribution inferred from 
Fe Il(A4570) /Fe Il(UV) 

Now we can roughly estimate column densities from the ob- 
served FeIl(A4570)/FeIl(UV) by using equation (5). The 
estimated column density distribution is shown in Figure 
8. An average and a standard deviation of the distribu- 
tion are found to be (log Nh , <J log n h ) = (22.8,0.5). Mar- 
coni et al. (2009) has suggested that radiation pressure does 
play an important role in BLR gas dynamics if column 
densities of BLR clouds have intrinsic dispersion such as 
(log Nh , criog n h ) = (23.0,0.5). Our results support that the 
assumption adopted in Marconi et al. (2009) is appropriate 
and that the radiation pressure plays an important role in 
BLR clouds. 



4.3 Correlation between Eddington ratio and 
Fe Il(A4570) /Fe Il(UV) 

Figure 9 shows the relation between FeIl(A4570)/FeIl(UV) 
and the Eddington ratio. A positive correlation is seen. 
Linear regression analysis, using an IDL procedure "FI- 
TEXY.pro" (cf. Press et al. 1992), gives the relation as: 



lof 



Feii(A4570) 
Feii(UV) 



-0.71 + 0.31 log 



(16) 



The Spearman's rank correlation coefficient 3 for assessing 
the nonlinear correlation is r$ — 0.58. This means the 



3 Spearman's rank correlation coefficient is non-parametric mea- 
sure of correlation, that is, which assesses how well an arbitrary 
monotonic function could describe the relationship between two 
variables without making any other assumptions 
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Figure 9. Eddington ratio vs. Fe II(A4570)/Fe II(UV) flux ratio. 
Typical lcr error is indicated at the lower right corner. The dashed 
line is the "FITEXY.pro" fit in the form of logj/ = a + /31ogx. 
The values of a, ft and Spearman's rank correlation coefficient rg 
are shown in the figure. 
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Figure 10. Estimated column densities vs. i^ol/^BddO- Each 
line represents L^oi = L^ddrad- Thick solid: a = 0.6 (adopted 
in Marconi et al. 2008); dotted: a = 1; dashed: a = 0.1; dot- 
dashed: a = 0.01; dot-dot-dot-dashed: a = 0.001. Note that a = 



Lion/ L^oi- All these lines intersect at JVjj 
Lbol/LEddfi = 1> where o T N H = 1. 



1.5 X 10 24 cm -2 and 



probability of the null hypothesis that there is no cor- 
relation is less than 10~ 13 . Thus the correlation between 
Fen(A4570)/Fell(UV) and the Eddington ratio is real. This 
implies that the column density increases with the Edding- 
ton ratio, because FeIl(A4570)/Fe Il(UV) increases with the 
column density. 

As was recently suggested by Dong et al. (2009), un- 
der the condition where the BLR clouds are subject to 
the radiation pressure, low-column-density clouds would be 
blown away by relatively large radiation pressure at large 
Lboi/LEdd.o, so that only high-column-density clouds would 
be able to be gravitationally bound. Figure 9 is a supportive 
evidence for their suggestion. 



samples becomes gravitationally bound. This means that 
the fraction of ionizing photons irradiating on Fe II-emitting 
clouds is much less than those on usual BLR clouds. Then 
it seems quite natural to conclude that the Fell emission 
does not originate in the region where usual emission lines 
such as H/3 originate, but originate in outer parts of BLR 
where the incident ionizing photon fraction becomes as low 
as a = 0.01. It is worth noting that from the studies of Ol 
and Ca II emission lines, M atsuoka et al. (2008) also suggests 
that Fe II emission originates in outer parts of BLR. 



4.4 Super-Eddington problem 

Figure 10 plots our samples on Nh — Lboi/ Lsdd,o plane. Each 
line represents L 00 i = LEdd,rad, so that the lower region of 
the line corresponds super-Eddington area. If we adopt ion- 
izing photon fraction a = 0.6 (i.e., thick solid line in Figure 
10), which is an average value for AGNs calculated by Mar- 
coni et al. (2008), almost all of our samples become super- 
Eddington. This result can be interpreted in two ways: (i) 
the conversion from Fe Il(A4570)/Fe Il(UV) to column den- 
sities (i.e, equation (5)) is wrong, or (ii) the adopted value 
of a is inappropriate. 

In the case (i), since FeIl(A4570)/FeIl(UV) is a func- 
tion of the column density and the temperature, the false 
in the conversion is attributed to the assumed temperature. 
If we assume hot clouds such as T e > 10, 000 K, the cor- 
responding column densities would become large, resulting 
in the solution for this super-Eddington problem. However, 
as already discussed, the previous studies favor cold clouds 
such as 6, 000 < T e < 10, 000 K for Fell emission (cf. Collin 
et al. 1980, Joly 1987). Thus this interpretation seems to be 
inappropriate. 

In the case (ii), as can be seen from Figure 10, if we 
adopt small values for a such as 0.01, the majority of our 



4.5 Eigenvector 1 in terms of the column density 

Boroson & Green (1992) applied a principal component 
analysis to low-redshift quasars and found that the prin- 
cipal component 1 (which is called "Eigenvector 1") links 
the strength of optical Fell emission and the weakness of 
[O III] emission. After a while, Boroson (2002) showed that 
Eigenvector 1 is driven predominantly by an Eddington ra- 
tio. However, the physical causes making up the Eigenvector 
1 has been left unknown. 

Here we propose a physical interpretation of Eigenvec- 
tor 1 in terms of the column density. As discussed in the pre- 
vious section, small- column- density clouds would be driven 
away from the line-emitting region by the radiation pres- 
sure at large Eddington ratio, and only large-column-density 
clouds can be gravitationally bound. Radiative transfer ef- 
fects make the optical Fell emission become large in such 
large-column-density clouds. On the other hand, ionizing 
photons emitted from the central object are intervened by 
these large-column-density clouds and thus have less proba- 
bilities of ionizing photons reaching to NLR clouds, resulting 
in weak [O III] emission. In fact, Figure 10(t) in Tsuzuki et al. 
(2006) shows a negative correlation between the [O III] /H/3 
and FeIl(01)/FeIl(Ul), which is almost the same as our 
FeIl(A4570)/Fen(UV), for 14 quasars. 
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5 SUMMARY 

(i) Analysis of the Fell(UV), FeIl(A4570) and Mgll emis- 
sion lines is performed for 884 SDSS quasars in a redshift 
range of 0.727 < z < 0.804. 

(ii) We suggest that FeIl(A4570)/FeIl(UV) can be an in- 
dicator of the column density of Fe II-emitting clouds regard- 
less of the excitation mechanism, i.e., photoionized or col- 
lisionally ionized clouds. From model calculations, we have 
confirmed this suggestion. 

(iii) Our photoionization models underpredict 
FeIl(A4570)/FeIl(UV) by a factor of 10, consistent 
with the preceding studies. Unless shocks are a viable 
heating mechanism, the failure of the photoionization 
model simply indicates that it is not predicting the correct 
heating rate, or that the radiative transport calculations 
are not correct. Ignoring anisotropy of UV Fell emission 
may be one of the causes. 

(iv) The column density distribution estimated from 
Fe Il(A4570)/Fe Il(UV) is almost the same as the one sug- 
gested by Marconi et al. (2009), supporting that the radia- 
tion pressure does work on Fe II-emitting clouds. 

(v) We also find a positive correlation between 
FeIl(A4570)/FeIl(UV) and the Eddington ratio, implying 
the links between the column density and the Eddington 
ratio. 

(vi) We find that under the assumption of the ionization 
fraction a — 0.6, almost all of our samples become super- 
Eddington. This problem can be cleared if the Fe II emission 
originates in outer parts of BLR where the ionizing photon 
fraction becomes as low as a — 0.01. 

(vii) We propose physical interpretation of 'Eigenvector 
1' in terms of the column density. In the interpretation, the 
strength of the optical Fe II emission results in the radiative 
transfer effects, while the weakness of the [Olll] emission 
results in the reduction of ionizing photons in NLR caused 
by intervening large column density BLR clouds. 
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